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Abstract. The study of detached eclipsing binaries is one of the most powerful ways to in- 
vestigate the properties of individual stars and stellar systems. We present preliminary masses, 
radii and effective temperatures for the eclipsing binary WWAurigae, which is composed of 
two metallic-lined A-type stars. We also reanalyse the data on HD 23642, an A-type eclipsing 
binary member of the Pleiades open cluster with a metallic- lined component, and determine 
its distance to be 139 ± 4pc. This is in agreement with the traditional Pleiades distance, but 
in disagreement with distance to the Pleiades, and to HD 23642 itself, derived from Hipparcos 
trigonometrical parallaxes. 
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1. Introduction 

The analysis of photometric and spectroscopic observations of double-Uned detached 
ecUpsing binaries remains one of the most powerful ways to investigate the properties 
of individual stars. The analysis of light curves and radial velocity curves allows the 
masses, radii and surface gravities of two stars to be found to an accuracy of better than 
1%. Further analysis using spectral synthesis or photometric index calibration techniques 
also allows the derivation of accurate effective temperatures, absolute magnitudes and 
distances of the components of the eclipsing system (Andersen, 1991; Southworth, Maxted 
& Smalley, 2004c). 

Accurate dimensions of a detached eclipsing binary (dEB) provide a discriminating 
test of theoretical stellar evolutionary models because the two components of the system 
have the same age and chemical composition but, in general, different masses and radii 
(Southworth, Maxted & Smalley, 2004b). The observation and analysis of dEBs in stellar 
clusters and associations may allow the calculation of accurate dimensions of more than 
two stars of the same age and initial chemical composition (Southworth, Maxted & 
Smalley, 2004a) . The distance of the cluster can be found from analysis of an eclipsing 
member, without the use of theoretical calculations and free of the difficulties of the main- 
sequence fitting technique. Furthermore, the age, metallicity and helium abundance of 
the cluster as a whole can be found from comparison of the observed masses and radii of 
the dEB with theoretical stellar models (Southworth et ai, 2004a). 

Another use of dEBs is to investigate the physical processes at work in single stars. 
This is particularly useful if one or both components belongs to a class of peculiar or 
poorly understood stars, for example slowly pulsating B stars (e.g., V539 Arae; Clausen, 
1996), or metallic-line stars (e.g., WW Camelopardalis; Lacy et al, 2002). Metallic A stars 
are well represented in the compilation of accurate eclipsing binary data by Andersen 
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(1991) but the physical processes and particular conditions of occurrence are still not 
fully understood. 

1.1. WW Aurigae 

WWAurigae (HD 46052, HIP 31173, V = 5.9, period = 2.52 d) is a dEB with an accurate 

Hipparcos parallax distance of 84±8pc (Ferryman, et al, 1997). Its eclipsing nature was 
discovered independently by Solvicv (1918) and Schwab (1918). Huffer & Kopal (1951) 
and Piotrowski & Serkowski (1956) undertook photoelectric observations but were both 
hampered by bright skies. Etzcl (1975) observed excellent photoelectric light curves in 
the Stromgren uvby filters, consisting of about one thousand observations in each filter. 

Kiyokawa & Kitamura (1975) published excellent UBV light curves of WW Aur and 
analysed them using a procedure based on rectification (Russell & Merrill, 1952; Kita- 
mura, 1967). Kitamura, Kim & Kiyokawa (1976) published good photographic spectra 
and derived accurate absolute dimensions of WW Aur by combining their results with 
those of Kiyokawa & Kitamura. 

The rotational velocities of the components were found to be 35 and 55kms~^ from 
CCD spectra by Abt & Morrell (1995), who also classified the stars as Am (A2,A5,A7) 
where the bracketed spectral types have been obtained using the Call K line, Balmer 
lines, and metallic lines. 

1.2. HD 23642 

HD 23642 (HIP 17704, V = 6.8, period = 2.46 d) was discovered to be a double-lined 
spectroscopic binary by Pearce (1957) and Abt (1958). Its Hipparcos parallax gives a 
distance of lll±12pc (Perryman, et al, 1998). Abt & Lcvato (1978) provided a spectral 
classification of AO Vp (Si) -|- Am, where the metallic-lino character of the secondary star 
relies on the presence of strong Fei lines. Torres (2003) discovered shallow secondary 
eclipses in the Hipparcos photometric data of HD 23642 and also presented an accurate 
spectroscopic orbit. 

Munari et al. (2004, hereafter M04) derived precise absolute masses and radii of both 
components from five high-resolution echelle spectra, and complete BV photoelectric 
light curves. From consideration of the effective temperatures and bolometric magnitudes 
of the components, and using the Wilson-Devinney light curve fitting code (Wilson & 
Devinney, 1971; Wilson, 1993) they found the distance to HD 23642 to be 131.9 ± 2. Ipc, 
where the quoted error is the formal error of the fit. 

HD 23642 is is an important system because it is a member of the Pleiades open cluster, 
for which the distance is controversial. 

1.3. The distance to the Pleiades 

The Pleiades is a nearby young open star cluster which has been exhaustively studied by 
many researchers. The 'long' distance scale of 132 ± 3 pc is supported by main sequence 
fitting analyses and photometric calibrations (e.g., Johnson, 1957; Crawford & Perry, 
1976; Meynet, Mermilliod & Maeder, 1993) and analysis of the astrometric binary Atlas 
(Pan, Shao & Kulkarni, 2004). Trigonometric parallax observations from the ground 
(131 ± 7pc; Gatewood et al, 2000) and from the Hubble Space Telescope (135 ± 3pc; 
Benedict, 2004) also agree with the 'long' distance scale. 

The 'short' distance scale of 120 ± 3pc is derived from Hipparcos parallaxes (van 
Leeuwen, 2004) but some researchers have suggested that the Hipparcos parallaxes are 
correlated on small angular scales (Makarov, 2002). However, vanLeeuwen (1999) has 
found other nearby open clusters with main sequences as faint as the Pleiades (located 
with the use of Hipparcos parallaxes). 
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Castellani et al. (2002) have shown that current theoretical steUar evolutionary models 
can fit the Pleiades main sequence if a low metal abundance of Z = 0.012 is adopted. 
However, Stello & Nissen (2001) used a metallicity- insensitive photometric technique to 
demonstrate that, if the Hipparcos parallaxes were correct, the main sequence Pleiades 
stars were implausibly fainter than their counterparts in the field. Also, Boesgaard & 
Friel (1990) have measured the iron abundance of the Pleiades to be approximately solar 
([§] ~ — 0.034 ± 0.024) from high-resolution spectra of twelve F dwarfs in the cluster. 

Pinsonneault et al. (1998) and Soderblom et al. (1998) investigated this discrepancy in 
distance determinations by performing a main sequence fitting analysis of the Pleiades 
and by attempting to discover nearby subluminous stars with solar metal abundance. 
They concluded that the Hipparcos measurements were most likely wrong and also noted 
that the most discrepant measurements belong to stars in the centre of the Pleiades, 
where the stellar density is the highest. This is consistent with the suggestion that the 
Hipparcos parallaxes are biased where the apparent stellar density is high. Further details 
can be found in Southworth, Maxted & Smalley (2004c). 

2. WW Aurigae 

2.1. Spectroscopic analysis 

A total of 59 grating spectra were obtained in 2002 October using the 2.5 m Isaac Newton 
Telescope on La Palma, covering the wavelengths 4220-4500 A at a two-pixel resolution 
of 0.2 A. 

Radial velocities were determined from the observed spectra using the two-dimensional 
cross-correlation algorithm todcor (Zucker & Mazeh, 1994). In this algorithm two tem- 
plate spectra are simultaneously cross-correlated against each observed spectrum for a 
range of possible velocities for each template. Template spectra of seven standard stars, 
with spectral types between AO IV and F5 V, were used, with circular spectroscopic orbits 
being fitting to the radial velocities calculated using each combination of spectra. The 
final spectroscopic orbital parameters are calculated from the mean and 1 a uncertainties 
of the individual good spectroscopic orbits. 

The final spectroscopic orbital parameters are Kp^ — 116.81 ± 0.23kms^^, Kb = 
126.49 ± 0.28kms~^, for a circular orbit with an orbital period of 2.52d. This gives 
a mass ratio oi q — 0.9235 ± 0.0027 and individual minimum masses of Ma sin'^ i = 
1.9590±0.0074 Mq and Mb sin^ i = 1.8090±0.0065 Mq. The best individual spectroscopic 
orbit has been selected and plotted in Figure ^ 

2.2. Light curve analysis 

The light curves of Etzel (1975), Kiyokawa & Kitamura (1967) and Huflfer & Kopal (1951) 
of WWAur were analysed using ebop (Nelson & Davis, 1972; Popper & Etzel, 1981), 
a simple and efficient light curve fitting code in which the discs of stars are modelled 
using biaxial ellipsoids. It is applicable only to eclipsing binaries which are sufficiently 
well detached for the stellar shapes to be close to spherical, as in the case of WW Aur. 
Investigations showed that the light curves display negligible eccentricity and third light 
so these quantities were fixed at zero for the final solutions. The light curves and their 
best fits are plotted in Figured 

Uncertainties in the values of the photometric parameters have been estimated using 
Monte Carlo simulations (Southworth, Maxted & Smalley, 2004b). After finding the best 
fit, the model light curve was evaluated at the phases of the observed datapoints. Random 
Gaussian noise, of the same size as the standard error of the fit, was added to simulate 
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Figure 1. Spectroscopic orbit for WWAur. Tlie observed radial velocities are shown with 
filled circles and the best-fitting spectroscopic orbit is shown using unbroken lines. 



WWAur A WWAur B 



Mass (Mq) 


1.964 


±0.007 


1.814 


±0.007 


Radius ( R©) 


1.980 


±0.009 


1.807 


±0.009 


Surface gravity logg (cms~^) 


4.160 


±0.007 


4.165 


±0.007 


Effective temperature (K) 


8350 


±200 


8170 


±300 


Equatorial rotational velocity (kms~^) 


35 


±10 


55 


±10 


Synchronous rotational velocity (kms^^) 


39.69 


±0.18 


36.24 


±0.18 



Table 1. Absolute dimensions and astrophysical parameters for the components of WWAur. 

observational noise, ebop was used to fit the resulting 'observation-like' light curve, and 
this process was undertaken ten thousand times for each real light curve. 

The uvby photometric indices for the individual components of WWAur are (b — y) = 
0.074±0.020and0.091±0.027,mi = 0.216±0.036 and 0.254±0.047 and ci = 0.977±0.039 
and 0.895 ± 0.0524 for the primary and secondary stars respectively. The calibration 
of Moon & Dworetsky (1985) then gives effective temperatures of 8350 ± 200 K and 
8170 ± 300 K for the two stars and indicates that the light of the system suffers from 
negligible interstellar absorption. If the system value of the H/3 index is assumed to be 
a reasonable estimate for the individual values for each star, the same results are found 
and it is confirmed that interstellar absorption is negligible. 

2.3. Absolute dimensions and comparison with stellar models 

The absolute dimensions of the two stars have been calculated from the results of the 
photometric and spectroscopic analyses and are given in Tabled The Granada theoretical 
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Figure 2. The Etzel (1975) difFerential light curves ol WW Aur, compared to the best-fittmg 
hght curves found using ebop. The residuals of the fits are plotted with magnitude offsets for 
clarity 

stellar evolutionary models (Claret, 1995) are able to match the astrophysical parameters 
of the components of WW Aur, for a solar chemical composition, for an age of 565 ± 
15Myr. The available predictions for other chemical compositions (Claret & Gimenez, 
1995; 1998; Claret, 1995; 1997) are unable to simultaneously fit the masses and radii of 
the two stars. Further details can be found in Southworth et al. (2004d). 



3. HD 23642 

3.1. Spectroscopic analysis 
A spectroscopic orbit was fitted to the five M04 echelle radial velocities for each star. The 



parameters of the orbit are 

,3 



Kb = 140.20 ± 0.57km s-\ giving 



99.10±0.58kms 

q = 0.7068 ± 0.0052, MAsin^i = 2.047 ± 0.016 Mq and A/b sin^ i = 1.447 ± 0.013 Mq. 
The orbital solution, plotted in Figure 01 is in reasonable agreement with those of Torres 
(2003), Abt (1958) and M04. 

From comparison between the observed spectra of M04 and synthetic spectra calcu- 
lated using UCLSYN and ATLAS9 model atmospheres (see Southworth et al, 2004a for 
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Figure 3. Comparison between the astrophysical parameters of the components of WW Aur 
and the Granada theoretical stellar evolutionary models (Claret, 1995). The predictions for 

metal abundances of Z = 0.01, 0.02 and 0.03 arc plotted using squares, circles and triangles, 
respectively. The isochrones have been plotted by interpolation in mass for all other quantities. 
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Figure 4. Circular spectroscopic orbit for HD 23642 derived using the Elodie radial velocities. 
The observed radial velocities are shown with filled circles and the best-fitting spectroscopic 
orbit is shown using unbroken lines. 

references), the effective temperatures of the two stars were found to be 9750 ± 250 K 
and 7600 ± 400 K, in good agreement with the results of M04. The uncertainties in these 
values include possible systematic errors caused by the slight spectral peculiarity of both 
components of HD 23642. 

3.2. Light curve analysis 

The B and V light curves of M04 were solved individually using EBOP. Although these 
light curves are of reasonable quality, deriving accurate parameters from them is problem- 
atic due to the shallow eclipses. As contaminating 'third' light, L3, is poorly constrained 
by the observations, we have made separate solutions for L3 = 0.0 and 0.05 (in units 
of the total light of the eclipsing stars) and included differences in the parameter val- 
ues derived in the uncertainties quote below. As there are no features in the spectra of 
HD 23642 known to come from a third star, it is unlikely that L3 is greater than 5%. 

The ratio of the radii, k, is also poorly constrained by the light curves; values between 
0.8 and 1.1 fit the observations equally well. To avoid this problem we have chosen to 
use the spectroscopic light ratio of Torres (2003), |j = 0.31 ± 0.03, found at 5187 A. 
This has been converted to a V filter light ratio using synthetic spectra, calculated from 
ATLAS9 model atmospheres, for the effective temperatures and surface gravities found 
in preliminary analyses. The resulting V filter light ratio of 0.335 ± 0.035 (where the 
uncertainties include a minor contribution due to possible systematic errors from the use 
of ATLAsQ model atmospheres) has been used to constrain k using the V light curve. The 
resulting values of k were then adopted for solution of the B light curve. 
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Figure 5. The M04 B and V light curves with our best fitted overplotted. The V light curve is 
shifted by +0.1 mag for clarity. The residuals of the fit are offset by +0.25 mag and +0.30 mag 
for the B and V light curves respectively. Note that the poor fit around the secondary eclipse 
in the B light curve is due to scattered data. 





Primary star 


Secondary star 


Mass (M0) 


2.193 + 0.017 


1.550 + 0.014 


Radius ( R©) 


1.831 + 0.030 


1.548 + 0.045 


Surface gravity logg (cms~^) 


4.254 + 0.018 


4.249 + 0.029 


Effective temperature (K) 


9750 + 250 


7600 + 400 


Equatorial rotational velocity (kms~^) 


38+1 


32 + 2 


Synchronous velocity (kms^^) 


27.7 + 0.6 


31.8 + 0.9 



Table 2. Absolute dimensions and astrophysical parameters for the components of HD 23642. 



The uncertainties in the fitted parameters were again estimated using Monte Carlo 
simulations, including significant perturbation of the linear limb darkeing coefficients for 
both stars. The best fit is plotted in Figure [S] 

3.3. Absolute dimensions and comparison with stellar models 

The absolute dimensions of the two stars have been calculated from the results of the 
photometric and spectroscopic analyses and are given in Table 12 Our resuhs are m 
reasonable agreement with the results of M04, and although the radius of the secondary 
star is somewhat larger, the two values are consistent within their uncertainties. 

In Figure ini the observed properties of HD 23642 have been compared with predictions 
of the Granada stellar evolutionary models for Z = 0.01, 0.02 and 0.03, adopting an age 
of 125 Myr for the Pleiades (Stauffer, Schultz & Kirkpatrick, 1998). The positions of the 
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Figure 6. Comparison between tlie observed properties of HD 23642 and the Granada models 
for metal abundances of Z = 0.01 (circles), Z = 0.02 (triangles) and Z = 0.03 (squares). 
Predictions for normal helium abundances are plotted with dashed lines and helium-rich model 
predictions (F = 0.34 and 0.38 for Z = 0.01 and 0.02 respectively) are plotted using dotted 
lines. An age of 125 Myr was assumed. 



components of HD 23642 in the mass-radius plane are consistent with the metal abun- 
dance of the Pleiades being solar or slightly supersolar. This confirms the atmospheric 
solar iron abundance determined by Bocsgaard & Fricl (1990) from high- resolution spec- 
troscopy of F dwarfs, but is in disagreement with suggestions that the 'short' and 'long' 
Pleiades distances could be reconciled by adopting a low metal abundance for the cluster. 

The components of HD 23642 cannot be fitted in the mass-radius plane by models of 
metal abundance or helium abundance significantly different from the solar abundances; 
for a solar chemical composition the predictions of the Granada stellar evolutionary 
models are consistent with the dimensions of HD 23642 for ages between 125 and 175 Myr. 
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Figure 7. The mass-radius diagram for A-type dEB components witii metallic-lined charac- 
teristics (filled circles) compared to other A-type dEB stars(open circles). The component stars 
of WW Aur are connected by a dashed line and the components of HD 23642 are connected by 
a dotted line. 

4. The Am star phenomenon 

There are a significant number of dEBs with components between 1.5 and 2.5 M0 
which have accurate masses and radii. We have taken a hst from Andersen (1991) and 
added eight systems which have been studied more recently (EI Cep, WW Cam, BP Vul, 
GGOri, FSMon, V364Lac and WW Aur and HD 23642). These stars have been plotted 
in mass-radius, temperature-gravity, period-mass and period-radius plots (Figures |H1 
ElandHnj. 

It can be seen from the figures above that there is no region of astrophysical param- 
eter space for which all stars are metallic-lined, although there is a cluster of Am stars 
around mass M = 1.9 and surface gravity logg = 4.15 (cms~^) which includes both 
components of WW Aur. Suggestions that metallic-lined stars have slightly larger radii 
than expected for their spectral type (Budaj, 1996) are therefore not confirmed. A-type 
stars of higher mass, though, do seem to need a lower surface gravity to be Am stars. 

As the predictions of theoretical models of approximately solar chemical composi- 
tion can successfully match the masses and radii of both HD 23642 and WW Aur, the 
metallic-lined spectrum phenomenon is confirmed to be a surface feature which does not 
significantly modify the interior conditions of the star. It is also clearly not a reflection 
of an overall high metal abundance. It remains unclear why some stars display metallic 
spectral lines when many similar ones do not. 

It should be remembered that using well-studied dEBs to investigate phenomena, such 
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Figure 8. The logarithmic effective temperature-surface gravity diagram for the sample of 
A- type dEB stars displayed in Figure Q Symbols have the same meaning as in Figure Q 



as enhanced metallic lines, can be misleading. The observational sample is very biased 
towards slightly evolved stars (Andersen, 1991) because these dEBs are more likely to 
show eclipses. A bias also exists towards dEBs with mass ratios close to unity because 
the observational requirements of such systems are more easily satisfied. 



5. The distance to the Pleiades 

5.1. Distance using bolometric corrections 

Adopting the astrophysical parameters of HD 23642 from Table El and using bolometric 
corrections from Bessell, Castelli & Plez (1998) to convert absolte bolometric magnitudes 
into absolute visual magnitudes, we find a distance to HD 23642 of 138.1 ±4.7pc. Whilst 
this is an acceptable result to adopt for the final distance to HD 23642, it is still affected 
by (probably small) systematic error due to the use of bolometric corrections calculated 
from theoretical models. 

The distance found by M04, 131.9 ± 1.9 pc, was derived using the above method but 
is significantly different from the result we obtain. The probable reason for this, which 
is consistent with the results given by M04, is that the absolute bolometric magnitudes 
given by M04 were calculated using values of the solar absolute bolometric magnitude and 
luminosity of Mboi© = 4.77 and Lq — 3.906x 10^^ W, whereas the bolometric corrections 
used (Bessell et ai, 1998) were calculated with Afboi© = 4.74 and Lq = 3.855 x 10^^ W. 
This inconsistency is sufficient to explain the difference between the distance derived by 
M04 and by ourselves, from the same observational data. 
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Figure 9. The period-mass diagram for the sample of A-type dEB stars displayed in 
Figure Q Symbols have the same meaning as in Figure Q 

5.2. Distance using surface brightness calibrations 

An alternative to considering the unobserved parts of the stellar spectral energy distri- 
butions is to use empirical relations between effective temperature and the wavelength- 
dependent surface brightness of a star. In this method the angular diameter of the star 
is compared to its linear diameter, determined from photometric analysis, to find the 
distance. The main source of uncertainty in this method is usually the light ratio, found 
from the photometric analysis, needed to calculate the individual magnitudes of the 
components of the dEB from the apparent magnitude of the system. If the effective tem- 
peratures of the two stars are known independently, this difficulty can be bypassed by 
comparing the apparent magnitude of the system to the flux-weighted combined angular 
diameter of the two stars. 

From consideration of the definitions of surface brightness and angular diameter, it 
can be shown that the distance to a dEB is given by 



A 



,(mA=0) 

''a 



+ 



2i?F 



,(mA=0) 



(5.1) 



where, for a distance given in parsecs, the stellar absolute radii, Ra and Rb, are given in 



,(mA=0) 



and 



are given in arcseconds 



,(m>>=0) 



AU, the zeroth-magnitude angular diameters 

and A represents a broad-band filter passband. Calibrations for 
terms of effective temperature by Kervella et al. (2004). The scatter of the fits is smaller 
at infrared wavelengths. We have applied the calibrations for the B and V filters (using 
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Figure 10. The period-mass diagram for the sample of A-type dEB stars displayed in 
Figure Q Symbols have the same meaning as in Figure Q 

Tycho apparent magnitudes transformed to the Johnson system) and for the JHK filters 
(using 2MASS data transformed to the SAAO system). We adopt the result from the K 
filter calibration, a distance of 139.1 ± 3.6 pc, as our final distance to HD 23642 and so 
to the Pleiades. The main source of uncertainty in this distance is from the uncertainty 
in the apparent K magnitude of HD 23642. 



6. Conclusion 

Preliminary results have been given from analysis of the A-type dEB WWAurigae, 
which is composed of two metallic-lined stars. The masses have been derived to accuracies 
of 0.3% and the radii to accuracies of 0.5%; predictions of the Granada theoretical stellar 
evolutionary models agree with the dimensions of the stars for an approximately solar 
chemical composition and an age of 565 ± 15Myr. However, the agreement between the 
uvhy and UBV light curves and the light curve of Huffer & Kopal is not good. If the latter 
data are rejected, the ratio of the radii increases to 0.953 ± 0.011, and the masses and 
radii of the two stars are unable to be fitted by theoretical models for one single age. For 
this reason, the study of WWAurigae presented here must be regarded as preliminary. 

HD 23642 is an A-type dEB which is a member of the Pleiades open cluster, for which 
the distance is currently controversial. The data of M04 have been reanalysed to provide 
robust parameter values and uncertainties. The metal abundance of the Pleiades has been 
estimated by comparing the masses and radii of the component stars of HD 23642 with the 
predictions of theoretical stellar models. The resulting metal abundance, approximately 
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solar, disagrees with suggestions that the 'short' and 'long' distances to the Pleiades 
can be reconciled by adopting a low metal abundance for the cluster. The distance to 
HD 23642 has been estimated to be 139±4pc from the use of surface brightness-effective 
temperature relations provided by Kervella et al. (2004). This distance is in agreement 
with the 'long' distance to the Pleiades and in disagreement with the Hipparcos distance 
to HD 23642 and the Pleiades. 
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